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ABSTRACT 


Analytical and experimental studies were initiated to determine 
if the response of a constant temperature hot wire anemometer to 
acoustic oscillations could serve as an analog to the response of the 
drop vaporization burning rate process to acoustic oscillations, and, 
perhaps, also as an analog to any Reynolds number dependent process. 

The motivation behind this study was a recent analytical study which 
showed that distorted acoustic oscillations could amplify the open-loop 
response of vaporization limited combustion. This type of amplification 
may be the cause of unstable combustion in liquid propellant rocket 
engines. 

The analytical results obtained to date for the constant 
temperature anemometer are similar in nature to those previously 
obtained for vaporization limited combustion and indicate that the 
response is dependent on the amount and type of distortion as well as 
other factors, such as sound pressure level, Mach number and hot wire 
tempe rature. 

In the experimental portion of the work, a system was designed 
and built which can create the same type of environment for the hot 
wire as was used in the theoretical analysis and which can determine 
the response for those conditions. Preliminary results indicate 
qua 1 i tat i ve agreement between theory and experiment. 



TABLE OF CONTENTS 


Page 

LIST OF FIGURES iv 

NOMENCLATURE vi 


Chapter 


1. INTRODUCTION 


1 


Setting 1 

Intent . . 2 

2. ANALYTICAL STUDY 4 


The Analog 

Hot-Wire Physics .... 
Flow Field Properties. . 
Nonlinear Response Factor 
Numerical Examples . . . 


4 

5 
11 
14 
16 


3. EXPERIMENTAL STUDY 


28 


Apparatus 28 

Procedures 40 

Results 48 


4. CONCLUDING REMARKS 


54 


REFERENCES. 


55 



LIST OF FIGURES 


Figure Page 

1. Schematic Diagram of the Hot Wire and Its Environment. ... 6 

2. Disa Hot-Wire Bridge Network 6 

3. Schematic Diagram of the Hot-Wire In Its Cylindrical 

Enclosure 12 

k. Nonlinear In-Phase Response Factor Versus Root-Mean- 

Square Sound Pressure Level for Various Harmonic 

Contents; M = 0.005, <J>i = 4>2 = 01*62= 0, 

T - T = 100 K° 19 

w °° 

5. Nonlinear In-Phase Response Factor Versus Harmonic 

Content for Various Sound Pressure Levels; 

M = 0.005, <h = <f>2 = 01 = 02 = °> T w " T «> = ,00K °- • • • 21 

6. Nonlinear In-Phase Response Factor Versus Root- 

Mean-Square Sound Pressure Level for Various 

Mach Numbers; d> 1 “ 4>2 33 e l ** e 2 ^ 0, T “ 100K°, 

P2/P1 - 0.7 W . • " 22 

7. Nonlinear In-Phase Response Factor Versus Mach Number 

for Various Harmonic Contents; 4>i = d>2 = e l = ®2 ~ 0, 

T - T = 100K°, P* = 0.00692 23 

w «. ’ rms 

8. Nonlinear In-Phase Response Factor Versus Hot Wire 

Overtemperature; M = 0.005, <j>i = cj > 2 = 01 = 0 2 “ 0, 

p 2 /Pi - 0 . 7 , P' = 0.00692 Zk 

ri rms 

9. Nonlinear Response Factors Versus Phase Angle Between 

Fundamental Sound Pressure Component and Its Second 
Harmonic; M = 0.005, 4>i = ©i = 62 "= 0, T - T = 100K°, 

P2/P1 = 0.7, P' - 0.00692 . W . . ” 25 

IlTO 

10. Nonlinear Response Factors Versus Velocity Phase Angle; 

M = 0.005, <h = $2 = 0, 0 2 = 0 i, T - T = 100K°, 

p 2 / Pl = 0.7, P' rms = 0.00692. . . . W . . " . 27 

11. Schematic Diagram of the Experimental System 29 

12. Photograph of the Experimental System - View 1 30 

13- Photograph of the Experimental System - View 2. 31 

iv 



V 


Figure Page 

14. Test Section With Right Angle Hot Wire Probe and 

Bare Microphone 34 

15. Sound Generation System 36 

16. Date Analysis System 38 

17- Hot Wire Calibration Equipment 41 

18. Hot Wire Calibration Curve 43 

19. Sound Pressure Variation Along a Diameter 45 

20. Change in Phase Along a Diameter 46 

21. Experimental and Computer Corsscorrelograms for the 

First Standing Tangential Mode 50 

22. Experimental Time-Delayed Crosscorrelation Ratio 52 

23. Predicted Time-Delayed Crosscorrelation Ratio for 

the Conditions of Figure 22 53 



NOMENCLATURE 


Sybmol 

A hot wire surface area 

w 

A, B hot wire constants 

A, B, C functions defined by (18) 

c speed of sound at T 

D diameter of test section 

D hot wire diameter 

w 

E hot wire bridge voltage 

E total energy 

h heat transfer coefficient 



hot wire current 

nonlinear out-of-phase response factor 
hot wire sensitivity to flow parallel to its axis 
thermal conductivity of air 
length of test section 
hot wi re length 
Nusselt number (hD w /k^) 
harmonic order (n = 1 ,2,3, - * -) 
air pressure 

undisturbed air pressure 

" 00 

rms pressure amplitude ■=• E p 2 

/ i n 
n=l 

harmonic pressure perturbation coefficient 
normalized cross correlation coefficient 



R gas constant for air 

d 

R , , cable resistance 

cbl 


vi 



Symbol 

Re 



u,v,w, 


V 


VV V z 


W 


*,y ,z 
e,r,z. 


Reynolds number (p VD /\i c ) 

1 g w f 

hot wire resistance at T 

o 

probe resistance 
standard resistance 
hot wire resistance at T 

w 

hot wire resistance at T 

00 

nonlinear in-phase response factor 
heat transfer 
air temperature 

film temperature, 1/2[T + T ] 

W \_ 

undisturbed air temperature 

harmonic temperature perturbation coefficient 
reference temperature 
hot wire temperature 
time 

cartesian components of velocity (x,y,z) 
wire cooling velocity 

cylindrical components of velocity (8,r,z) 

harmonic transverse air velocity coefficient 

transverse air velocity 

work transfer 

cartesian coordinates 

cylindrical coordinates 


Greek Symbols 

a hot wire resistance coefficient 

a hot wire angle of inclination to radius vector 

a 1Q eigen value (0.5861) 
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Greek Symbols 

Y 

G 

n 

P 

P„ 


♦ 


n 


U) 


ratio of specific heats 

harmonic phase angle for velocity-pressure relation 
ai r densi ty 

harmonic density perturbation coefficient 
harmonic phase angle for pressure perturbation 
air absolute viscosity 
radian frequency 


Superscripts 

1 time-dependent component 

time-mean component 
dimensionless property 


Subscripts 

f 

o 


evaluation at film temperature 

val ue at T = T = 0°C 
o 

undisturbed value 



Chapter 1 


INTRODUCTION 

Setting 

The need exists to obtain a more accurate understanding of the 
nonlinear instability limits of rocket engines. These limits are 
currently one of the least understood properties of combustion insta- 
bility. The causes of instability and the factors determining the 
instability limits are of major importance to the design of high 
performance rocket engines. 

Steady combustion in liquid propellant rockets is assumed to be 
vaporization limited. Priem and Heidmann [1] have shown that quantita 
tive agreement exists between theory and experiment. Vaporization 
limited combustion is also expected for unstable combustion and this 
has been qualitatively confirmed. 

Recently, Heidmann [2,3] has shown analytically that distorted 
acoustic oscillations affect the open-loop response of vaporization 
limited combust ion--a velocity sensitive, or Reynolds number dependent 
process. In particular, the open-loop response factor increased as 
the level of distortion increased. Response factors an order of 
magnitude greater than that obtained for sinusoidal excitation were 
found for some conditions. 

The vaporization model used in references 2 and 3 assumed that 
the burning rate was proportional to the drop Reynolds number raised 
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to a power. Such a model characterizes Reynolds number dependent 
processes other than drop vaporization burning rates. Thus, the 
results of references 2 and 3, and of this research, may also be of 
significance to convective heat and mass transfer processes in general. 

Intent 

The objective of this research is to experimentally confirm 
the results of Heidmann's studies by determining the open-loop response 
of a Reynolds number dependent process which is an analog of the drop- 
vaporization burning-rate process. An analogous process, which lends 
itself to direct measurement of the flow quantities, will greatly 
simplify the work over that of attempting to make measurements in a 
burning mixture. The ultimate objective, towards which this research 
is directed, is to determine the closed-loop response of the drop 
varpoization process to distorted acoustic waves; i.e., to determine 
if rocket combustion will become unstable when excited by distorted 
waves . 

The work to confirm Heidmann's results can be categorized as 
fol lows : 

•Engineering analysis of the open-loop response of the analog 
process to distorted waves. This is necessary to meaningfully 
interpret the experimental results. 

•Design and construction of a test system to measure the in-phase 
open-loop response factor for the analog process. 

• Measurement of the open-loop response factor for spinning trans- 
verse acoustic oscillations generated by pure sinusoidal and 
distorted sinusoidal input signals. This mode of oscillation 



would be studied since It is the most frequently encountered mode 
In unstable combustion [3]. Variations in the through-flow 
velocity, the harmonic amplitudes and the harmonic phase angles 
would be considered. 

• Comparison of the measured response factors with those predicted. 

At this time the first two items have been largely completed. 
The last two items are still in the preliminary stage. Each Item will 
be separately discussed. 



Chapter 2 


ANALYTICAL STUDY 
The Analog 

A qualitative description of the drop vaporization process will 
provide a basis for choosing a suitable analog. First we have a liquid 
droplet at a temperature close to its saturation temperature corre- 
sponding to the partial pressure of its vapor. The gas mixture, which 
surrounds the droplet, moves relative to it and has a temperature which 
is often significantly larger than the droplet's. Much of the energy 
transfer from the gas mixture to the droplet is by forced convection and, 
under steady-state conditions, results in vaporization of the droplet 
1 iquid. 

Assuming the drop vaporization process to be convection dominated, 
a reasonable analog of this process is that of forced convection heat 
transfer from a hot wire to a cool gas. Since the constant-temperature 
hot-wire anemometer provides an excellent means for experimentally 
determining the instantaneous energy transfer rate from a hot wire, it 
is being used as the analog process. 

To provide a controlled unsteady environment for the hot wire, a 
resonant sound field with forced mean gas motion will be used. Before 
giving a detailed description of the wire's environment, the physics of 
a constant-temperature hot wire in an unsteady gas flow will be developed 


k 
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Hot-Wire Physics 

Physical System 

The physical system is shown in figure 1. It is assumed there 
exists a steady mean flow, U, and that the wire is oriented such that 
its axis is normal to the mean flow direction. For the cartesian 
coordinate system shown in figure 1, the x-axis is in the mean flow 
direction and the z-axis is parallel to the cylinder's axis. 

Energy is electrically supplied to the wire such that its 
temperature is constant — the constant-temperature anemometer provides 
the needed energy and controls. A natural limitation for the anemometer, 
however, is that the gas temperature must not exceed the normal operating 
temperature of the wire. 

Energy Analysis 

Taking the hot wire to be a closed system, the first law of 
thermodynamics is 


dQ dW _ d£ 
dt " dt ~ dt * 

For constant temperature operation the total energy of the wire is 
constant and (I) becomes 


dQ _ dW 
dt “ dt * 


The rate of electrical work done by the system is 


dW 

dt 



( 2 ) 


( 3 ) 







7 


Three modes of energy transfer as heat exist. However, radiation 
is small enough to be neglected, conduction and convection are lumped 
together, and Newton's law of cooling is used to express the heat transfer 
rate; that is, 

4£«-hA[T -Tj. (b) 

at w w “ 

Introducing the Nusselt number, 

hD 

equations (2), (3) and (b) reduce to 

|2 R = rrk f L Nu f [T - T ] . (6) 

w r w r w °° 

Auxi 1 iary Relations 

Certain auxiliary relations are now introduced to permit (6) to 
be expressed in terms of the properties normally measured and/or known. 

The electrical resistance of the wire is practically a linear 
function of temperature. Expressed mathematically, this is 

R = R + a[T - T ], (7) 

wo 1 w o * 

where R is the wire resistance at T — T = 0°C in this analysis, 
o o o ' 

Letting R^ denote the wire resistance when it is in thermal equilibrium 
with its environment, the wire over-temperature is 



In the 01 SA anemometer system, the potential drop across the wire 
is not directly measured, and it is necessary to express the current in 



8 


(6) In terms of the DISA bridge properties. For the bridge network shown 
in figure 2 (on page 6), the hot-wire current is 


I = 


R + 
s 


cbl 


R + 
P 


(9) 


w 


There are a number of possible correlations that could be used 
to express Nu^. in terms of the flow field variables. One of the most 
reliable correlation forms, for low Reynolds number flows (0.02<Re^<4^) , 
seems to be that of the Col 1 is-Wi 1 1 i ams equation [4], namely 

Nu f = [A + B Re f 0 * 45 ] [T f /Tj 0 * 17 , (10) 


whe re : A , B 



dimensionless constants (A is a strong function 

of wire length when L /D < 400 and B seems to 

be a weak function of L /D .) , 

w w 

wire Reynolds number, V D w /Pf, 

wire cooling velocity, [(u + u') 2 + (v 1 ) 2 + K(w') 2 ] 0 ’ 5 , 
property evaluated at the film temperature, and 
film temperature, 1/2 [T + T ]. 

W 00 


Since the experiments will use atmospheric air at room temperature 
and pressure, the viscosity and thermal conductivity data will be taken 
from Hilsenrath [5] and the density will be computed from the ideal gas 
relation, namely, 


R 

a 


287.06 


m 


sec - K 


with 
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The Hilsenrath data in the temperature range of interest (273.15°K < T < 
**73. 15°K) , can be accurately approximated by 


and 


wi th 


ko (TA 0 ] 0 - 877 , 


y o lT/T o 10 ' 757 ’ 


02 ) 


(13) 


k 

o 


2.*»13 x 10" 2 


kg-m 

sec -K° * 


v = 1.716 x 10 
o 


-5 kg 


sec-m 


and 


T q = 273* 15°K. 


Bridge Voltage - Flow Property Relationship 

The relationship between the hot-wire bridge voltage and the flow 
properties is obtained by combining (6), (7), (9), (10), (11), (12) and 
(13). This gi ves 


E 2 - 


[R s + 


cbl 


R + 
_ E 


R ] 

w 


w 


[wL k 3 [/\ + B Re f 0.45(^2)0.225]. 

o 


0.17 


T 

i- o 


0.877 


[T - T ] , 
w “ 


where 


R = R + a[T - T ], 

W O W 0 


t- - 0.5 It + t 3 , 

f W a> 


m 


(7) 


( 15 ) 
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Re 


P V 0 
» o w 


f ° Vf M o tT f /T o ]0- 757 


V = 
o 


m 

sec 


V 2 = [ ( u + u 1 ) 2 + (v 1 ) 2 + K(w') 2 ]/V 2 , and 

o 

K - 0. 


(16) 


07 ) 


It was assumed that the wire is not cooled by fluid motion parallel to 
its axis and thus K was set equal to zero. 

Since the empirical relationship for Nu^, equation (10), was 
obtained for steady-state, steady-flow conditions, its applicability 
to unsteady flow is questionable. However, it is customary to assume 
that it holds for unsteady flows if instantaneous values are substituted 
for the fluid properties. This is done here and it constitutes a major 
assumption. 

Noting that the leading coefficient in (14) is constant during 
normal operation, and that it depends on the wire temperature, (14) can 
be rewritten in a simpler form by defining dimensional forms of A and B, 
namely, 


where 


A - C(T ) A, 
w 


B = C(T )B, 
w ’ 


(18a) 

(18b) 


C(T ) 
w 


trL k 
w o 


[R + R , , + R + r ]‘ 
1 s cbl p w 

R 

. w 


(18c) 
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With this nomenclature the final form of (14) becomes 




rTj 

0.17 

r T. 1 

0.877 

E = 

[A + B Re^* 45 ^ 2 ) 0 - 225 ] 

f 

T 


f 

T 

[T -TJ 
w 00 


. 

_ 00 


L o J 

— 


To determine E, the flow field properties must be known. 


0.5 

. (19) 


Flow Field Properties 

The flow field providing the wire's environment is schematically 
shown in figure 3- The cylindrical enclosure, with inside diameter D 
and length L, is acoustically closed but permits a uniform forced flow 
in the axial direction. 

With proper excitation of this enclosure, acoustic velocity 
oscillations in any of the three vector directions can be obtained. 
However, Heidmann [3] restricted his analysis to transverse-velocity 
oscillations — "because of the greater interest in transverse-mode 
instability in rocket combustors." Considering only transverse 
oscillations, the hot-wire cooling velocity can be expressed as 


V 2 = u 2 + (v') 2 . (20) 


Noting that V 2 is given in terms of velocity components based on the 
wire's cartesian coordinate system, it is necessary to relate these 
velocity components to those based on the enclosure's cylindrical system. 
The wire is oriented (figure 3) such that its axis is in the r-0 plane 
and makes the angle a with the radius vector. Thus 


u = v z > 


v'=v sin o + v, 


( 21 ) 

( 22 ) 


cos a 
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and (20) becomes 

wi th 
and 

All the other properties that are affected by the acoustic 
oscillations, thereby having acoustic components, will be expressed in 

t 

the form 

0 = 6 + 6 ' ( 24 ) 

* 

or 

6=8 [1 + 0 * 3 , ( 25 ) 

*■ 

where B is the time-mean component of the general property B, B‘ is its 
acoustic component, and 8' is its "normalized" acoustic component -- 
B' = B'/B. 

The expressions used for the normalized acoustic components are 
the same as those in reference 3» namely 

P' 

p ' = E p cos (no)t - <k ) , 

, n n 

n=l 

co 

T 1 = E T cos (nu>t - <f> ) , 

, n n 


00 

= E p cos (nu>t - d> ) . 
, r n y n 

n= 


_ 2 _ A A , 

V 2 =v [1+v sin a + v„ cos a r ] , 

z i r q 


( 23 ) 


V = v /v , 
r r z ’ 


u e ' v e / 5 z- 
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oo 

v = Z V cos (ntot - <J> - 0 ) , 

t . n n n 

n=l 

V t " v r ’ 

and 

v t = v t^z' ( 2 6) 

The expression for T' has been added to those given by Heidmann. 
The phase angle <j> specifies the phase relation between harmonic 

I 

components- It has been assumed to be identical for acoustic pressure, 

density and temperature- The gas velocity expression also contains <j> n 

but it is further modified by a velocity phase angle 0 n which specifies 

the velocity-pressure phase relations. When © n is zero, the velocity 

and pressure are in phase and traveling wave properties are simulated. 

A 0 of 90° simulates standing acoustic modes, 
n 

Nonlinear Response Factor 

Exposing a constant-temperature hot wire to a mean flow perturbed 
by acoustic oscillations results in an unsteady hot-wire energy transfer 
which is exhibited as an unsteady hot-wire voltage E. Expressed in terms 
of its mean and fluctuating components, E is 



E = E + E' , 

(27) 

or 




E = E [1 + E ■] , 

(28) 

with 




E' = E'/E. 

(29) 


A parameter which compares that portion of E' which is in phase 
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with the acoustic pressure P 1 is the nonlinear in-phase response factor 
defined by 


. 2u 


I E * (tot) P'(u)t) d(wt) 

3 = 2 

nl r 2n 

J [P 1 (cut )]2 d(ut) 

o 


(30) 


A parameter which compares that portion of E' which is in phase 
with the 90 -degree-phase-shifted acoustic pressure P' s is the nonlinear 
out-of-phase response factor defined by 


-2tt 


J E'(u)t) P* (wt) d (tot) 

_ o _ 

nl ~ rlv 

J [P’ s (u>t)l 2 d(uit) 


(3D 


where 


P' (wt) = E p cos(nwt - <(> - ir/2) , 


ps 


n=l 


or 


P' (wt) = E p„ sin(nwt - 4> ) . 


PS 


n= 1 


(32) 


Equations (30) and (31) are based on the approximation given by 


E'(u)t) - P'(a)t) + 1^ Pp S ( w t) . 


(33) 


To evaluate R nl and 1^, E'(o)t) must be related to the sound field 
given by equations (26). This is done in the following manner: 
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• The hot wire's properties, including i ts , operat ing temperature, 
are specified. 

•The hot wire's environment is specified. This includes the 
time-mean through-flow velocity v z> the time-mean absolute temperature 
and pressure of the air, and the sound field. 

•The instantaneous value of each property is concomitantly 
substituted in (19) » thereby producing E(ajt) . 

• The time-mean value of E(cat), E, is obtained from 


-2tt 


E = 


2^ I E(wt) d(a)t) 


(3*0 


Due to the complex mathematical form of E(u)t) , the integration is 
normally performed numerically. 

• The desired E'(u>t) is given by 

E ' (cut) = [E(u)t) - E]/E . (35) 

With E'(u)t) given by (35) and P'(u>t) given by (26), and 
are obtained from (30) and (31), respectively. The nature of the 
nonlinear response factors is best illustrated by numerical examples. 


Numerical Examples 

The numerical examples to be presented are based on the following 
hot-wire and environmental properties: 

• Wire and Bridge Properties 

L = 1 .02 mm 
w 

D = 5 pm 
w 
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a = 0.0129 ohm/K° 

R = 0.40 ohm 
P 

R k | =0.30 ohm 
R g = 100 ohm 

R w = [0.0129 T w (°K) - 0.08] ohm 

• Basic Air Properties 

T o =273.15°K 

k = 2.-413 x 10 -2 kg-m/(sec 3 -K°) 
o 

= 1.716 x 10 5 kg/(sec-m) 

R = 287.06 m 2 /(sec-K°) 

3 

y = 1.4 

V = 1 m/sec 
o 

ft 

• Wire Environmental Properties 

T, = T - 295.37°K 
P - P - 738.9 mm Hg 
L = 30.48 cm 
D = 20.42 cm 

With these properties specified, the main variables are: the 

hot-wire Mach number M = v /c , the hot-wire temperature T , and the 

z ® w 

acoustic properties. 

Traveling Acoustic Waves 

The examples for traveling waves will be for the first traveling 
transverse mode of resonance with the wire oriented such that 0=0 
and a = 0 (the wire's axis is parallel to the radius vector). For radial 
positions outside the acoustic boundary layer region, the acoustic 
density, temperature and particle velocity are simply related to the 
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acoustic pressure, assuming the harmonic components of the acoustic 
properties are related to those of the acoustic pressure in the same 
manner that they are for undistorted waves. With this assumption the 
harmonic components of equations (26) are related to the p n in the 
following manner: 


p n ~ Y 


and 



v 

n 


I R/r 

Y ttcx 10 P n’ 


6 = 0 . ( 36 ) 

n 

Heidmann's [3] velocity coefficient differs with that of (36) by its 
sign and the radially dependent coefficient R/(ira 10 r). The sign 
difference has no effect on the response factor. For purposes of 
comparison, the radial position will be taken as that for which R = naio r 
Thus the harmonic velocity components are given by 


v 

n 


1 

Y P r 


(37) 


Influence of Harmonic Content on R^. Considering sound pressure 
wave distortion due only to the second harmonic of pi (i.e., due to P 2 ) , 
the nonlinear in-phase response factor was computed. The results are 
displayed graphically in figure 4 for various ratios of p 2 to pj. For 
this Mach number (M = 0.005), maximum response seems to occur for a 
normalized root-mean-square sound pressure level of 0.00692 and for P 2 /P 1 



Figure 4. Nonlinear In-Phase Response Factor Versus Root 
Mean-Square Sound Pressure Level for Various 
Harmonic Contents; M = 0.005, <f>i = <t >2 = e l = 9 
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approximately equal to 0.7* Figure 4 is similar in form to Heidmann's 

[3] figure 2c. Note, however, that peak response occurs at a lower 

value of P 1 for the hot-wire, 
rms 

The value of P 2 /P 1 for which is a maximum (0.7) is clearly 
demonstrated in figure 5. From the three curves presented it is seen 
that this value is not dependent on the sound pressure level. 


Influence of Mach Number on R^. For a fixed pressure waveform 
(P 2 /P 1 = 0.7), R^ was computed for various wire Mach numbers. The 
results are shown in figure 6. Generally speaking, the lower the Mach 
number, .the greater the response. This is more clearly shown in 
figure 7, where three different harmonic contents are considered. 

It should be noted, in figure 6, that the value of the sound 
level at which peak response occurs increases with increasing Mach 
number. 

Comparing Heidmann's [ 3 ] figure 7 with figure 6 it Is seen that 
they agree in form. 

Influence of Wire Overtemperature on R^. The influence of 
wire overtemperature on the peak response for the conditions of figure 
4 is now considered. Fixing the sound field conditions (P 2 /P 1 = 0.7, 
P rms = 0 * 00692 ), the effect of wire overtemperature on R nl is shown 
in figure 8. Note that reaches its maximum value at about 
T w " T ® = i 75 K° . 


Influence of Harmonic Phase Angle on 1^ and R^>. The influence 
of <j >2 on the nonlinear response factors for the conditions of peak 
R nZ oP f'9 ure **« except for 4> 2 > is shown in figure 9* Note, that 
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M = 0.002 



P' 

rms 


Figure 6. Nonlinear In-Phase Response Factor Versus Root- 

Mean-Square Sound Pressure Level for Various Mach 
Numbers; <j>i = $2 * 8 l = e 2 = = 100 K° , 
p 2 /pi = 0.7 • 






2 


0.0k 


nl inear In-Phase 

r Various Harmoni 

- T - 100 K°, P 
00 7 
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is zero for <J» 2 = 0 and 180 degrees and that it takes on its maximum 
positive value at 4> 2 = 90 degrees. The in-phase response factor has 
its maximum at <j> 2 = 0 degrees. 

Figure 9 is similar to Heidmann's [3] figure 6. It is a cut 
of his figure at peak response. 


Influence of Velocity Phase Angle on 1^ and R^. By again' 
using the sound field properties that produced the peak value of R ^ 
in figure 4, the effect of 0 on the nonlinear response factors was 
determined and the results are shown in figure 10. 

The in-phase factor in this figure represents a cross- 
correlogram of E‘ and P‘ and the out-of-phase factor 1^ is a cross- 
correlogram of E' and P' s . Both of these can be obtained experimentally. 

Figure 10 represents a cut of Heidmann's figure 5 at peak 


response. 
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0 0.2tt 0.4n 0.6n 0.8n . 1 .On 


0 1 

Figure 10. Nonlinear Response Factors Versus Velocity Phase 
Angle; M = 0.005, $1 = $2 ~ 0> ®2 = ®l> 
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Chapter 3 


EXPERIMENTAL STUDY 

The object of the experimental portion of the research is to 
determine the open loop response of a constant temperature hot-wire 
anemometer to distorted acoustic oscillations, and to compare the 
results with those obtained theoretically. An experimental system is 
needed which will create the same velocity and acoustic environment 
as that used in the analytical portion of the work. It is also desirable 
to incorporate maximum flexibility into the system so that it can be 
readily adapted to any future work. 

The system needed consists of four basic parts — the test 
section itself, a means of generating the necessary acoustic oscillations, 
a means of analyzing the data, and the various types of equipment needed 
to monitor the performance of all parts of the system. Figure 11 is a 
schematic illustration of the system and figures 12 and 13 show two 
different views of the completed system. 

Apparatus 

Test Section 

The test section is an aluminum cylinder with a 20.42 cm I.D., 

25.4 cm O.D., and an effective inside length of 30.48 cm. The ends of 
the cylinder are closed with 1.27 cm thick circular aluminum plates into 
which approximately 310 equally spaced 2.54 mm diameter holes have been 
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Figure 11. Schematic Diagram of the Experimental System. 
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Figure 12. Photograph of the Experimental System - View 
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drilled so that a flow of air can pass through the section. The 
cylinder and the endplates have an 0-ring seal between them and the 
whole system is clamped together with four long bolts. The cylinder 
is divided into two sections, one 12.70 cm long, the other 17-78 cm 
long. They are joined with a lap joint containing an 0-ring seal. The 
short section is for the acoustic drivers; the long section is for the 
hot wire and microphone probes. By dividing the cylinder into two 
parts, the drivers and probes can be positioned at any angle relative 
to each other. 

The shorter cylindrical section has three threaded holes located 
at 0°, 90°, and 180° around the circumference, 6.35 cm from either end. 
Acoustic drivers can be mounted in any or all of these holes and plugs 
are available to close the holes when not in use. How many drivers 
are used, and which holes they are mounted In, depends on the type of 
sound field desi red. 

On the longer cylindrical section, provision, is made for mounting 
two probe holding positioners, one which is adapted to hold a 1.27 cm 
condenser microphone with a probe tube extension, and one which can hold 
the hot wire probe support. Both have micrometer type scales which 
can position a probe to within 0.25 mm. The microphone holder is located 
so that the microphone probe moves along a diameter at the center of the 
test section. This enables a radial pressure profile to be measured. 

When it is only necessary to know the pressure at the wall of the chamber, 
the microphone positioner can be removed and a bare 1.27 cm microphone 
can be flush mounted in the wall of the chamber. Using a microphone 
without a probe tube extension eliminates the complex, frequency dependent 
calibration curve associated with it. 
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The hot-wire positioner is designed to hold Disa miniature 
hot wire probe holders. Its center is located 7 cm from the center of 
the test section so that when a right-angle probe holder is being used, 
the hot wire is located at the center of the test section and can be 
moved along the same radial line as the microphone probe. The right- 
angle probe is relatively short and can only penatrate a distance of 
6.25 cm into the test section. To obtain velocity readings across the 
entire diameter of the section, a long straight probe holder can be 
used. However, for this case the hot wire will not be moving along a 
diameter at the center of the section. Figure 14 is a photograph of 
the test section with a right-angle hot wire probe holder in the system 
and a bare 1.27 cm condenser microphone flush mounted in the wall. 

To generate a flow of air through the test section, it is 
connected via a sheet metal cone and flexible tubing to a fan which 
draws air through the system. The cone, flexible tubing, and fan 
housing can be seen in figures 12 and 13. To keep from drawing dust 
into the system, which might alter the characteristics of the hot-wire, 
or destroy it, an air filter was mounted at the inlet to the test section. 
This had to be eliminated, however, since access to the inlet had to be 
maintained in order to check the sound field being generated. 

Sound Generation System 

The distorted acoustic oscillations, which are required, consist 
of a spinning tangential cosine wave at the fundamental frequency with 
a cosine wave at twice this frequency added to it. The intensity of 
this sound should be about 160 dB or greater so that the acoustic 
pressures and velocities are significant compared to the background noise 
in the system. 
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Test Section with Right-Angle Hot Wire Probe and Bare Microphone. 
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Theoretically, to create a spinning mode, two drivers are required, 
located 90° apart, which generate signals of the same amplitude and 
frequency, but 90° out of phase with each other. In actual practice, 
the signals going to the drivers may have to have slightly different 
amplitudes, and have a phase difference other than 90°, to produce a 
spinning mode. This is because the response of the drivers may not be 
exactly the same and the chamber itself interacts wi th the drivers. Thus, 
to generate the type of acoustic oscillations desired, two signals will 
have to be created containing the fundamental and the second harmonic 
frequencies, but whose amplitudes, both fundamental and second harmonic, 
can be adjusted independently, and whose phase differences can be 
adjusted relative to each other. Then, to create the desired sound 
field, adjustments to the phases and amplitudes can be made while 
monitoring the resulting sound field with an auxiliary microphone. The 
system used was outlined in Figure 11 and is shown in more detail in 
Figure 15. A cosine wave generated by the General Radio 1209"A function 
generator and tuned to the fundamental frequency is used as the basis of 
the sound generation system, with all other equipment adjusted relative 
to i t . 

The acoustic drivers are rated at 75 watts, but experience has 
shown that they cannot be used at this level when being driven by a 
pure tone, or a pure tone with one of its harmonics, without overheating. 
Because of this, they should not be driven to power levels above about 
50 watts each. However, this is sufficient to generate sound levels of 
over 160 dB. When less complex acoustic oscillations are desired, parts 
of the system are simply not used. 
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Data Analysis System 

The signals from the hot-wire and microphone must be processed 
to obtain a crosscorrelogram from which a response factor can be obtained. 
The system used was outlined in figure 11 and is shown in more detail in 
figure 16. 

The signal from the hot-wire is processed by the constant temper- 
ature anemometer. The output of the anemometer is the sum of two 
components — one caused by the through flow through the test section, 
the other a result of the acoustic velocity. An oscilloscope is used 
to remove the time-mean component of the anemometer signal. A sound 
level meter converts the output of the microphone to a signal proportional 
to the sound level. Both the resulting signals, E 1 and P', are then 
recorded on tape so the data is permanently recorded and can be analyzed 
from different aspects at later times. Also, recording the signals is 
in part made necessary by the intensity of the sound created by the 
system while in operation. It cannot be tolerated for long periods of 
time; hence, it is better to record the data, turn off the system, and 
then analyze the data. 

Phase shifts and/or time delays may occur in the signals during 
processing. Because of this it is necessary to incorporate a time delay 
unit into the system to compensate for them. The time delay unit is 
controlled by the level of the DC voltage input to it. Two inputs are 
available to the system. One, a DC power supply, can be adjusted 
manually to the particular voltage and, hence, time delay, that is 
desired. The other is a function generator which inputs a continuously 
varying DC signal which can be adjusted to scan the entire range of 
time delays from zero to a complete fundamental period. 
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Figure 16. Data Analysis System 
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After the signals have been delayed with respect to each other, 
they are each adjusted to 0.5 volts rms and added. A correlation coeffi- 
cient is then calculated from this sum. This correlation coefficient 
can be processed in two ways. If the time delay is a constant (DC power 
supply controlling the time delay), an integrating digital voltmeter 
and a timer can be used to obtain a time averaged correlation coefficient 
for that particular time delay; or, if the time delay is being 
continuously changed, the correlation coefficient can be plotted against 
this time delay on an x-y recorder to obtain a crosscorrelogram of the 
data. 

Auxiliary Monitoring Equipment 

To operate the system properly, and to be able to analyze the 
data, various auxiliary monitoring equipment is necessary. The only way 
to determine what type of sound field is being generated is to monitor 
it with a microphone. A second 1.27 cm condenser microphone with probe 
tube extension is part of the system. The probe tube extension is 1 mm 
in diameter and thus can be inserted into the system through the holes 
in the endplates. In this way the sound pressure level can be checked 
at various angular positions at a constant radial distance from the 
center; or, the pressure variation along a diameter can be checked. 

A phase angle meter is used for checking the phase between 
various signals. A dual trace oscilloscope is available to view the 
shapes of the various signals and especially to monitor the wave forms 
generated by the microphone and anemometer during data taking. To 
determine the frequency of the fundamental and its second harmonic, a 
frequency meter is available. These frequencies are not necessarily 
constant since they vary with temperature. Also, parts of the system 
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are various DC and AC voltmeters, wet and dry bulb thermometers, and 
a barometer. 


Procedures 

Before data is taken, all the equipment is adjusted and calibrated. 
Then the through flow and acoustic field are generated and data is 
recorded from the anemometer and microphone outputs. With the sound 

off, a crosscorrelogram is then made from the data. From this, response 

h 

factors can be calculated. Different types of acoustic fields were 
generated in order to better understand the system itself and to help 
in analyzing the data for the spinning transverse mode. 

Wi re Cal i brat ion 

To make a meaningful comparison between the analytical and 
experimental results, the equation relating air velocity to voltage 
output from the anemometer, which is used in the analysis, must be the 
equation that fits the particular wire being used in the experiment. 

To determine this equation, the wire has to be calibrated. The equip- 
ment used for this purpose is shown in figure 17 and consists of a very 
low turbulence wind tunnel and a micromanometer. These are used in 
conjunction with the hot wire anemometer to obtain the calibration. The 
tunnel has provisions for the mounting of a wire, connections to attach 
the wire to the anemometer and taps for connection to the manometer. 

The velocity of air through the tunnel can be determined from the pressure 
drop indicated on the manometer. To calibrate the wire, an overheat 
ratio must first be decided on. This determines the operating temperature 
of the wire relative to the ambient temperature and its value is a 




Figure 17- Hot Wire Calibration Equipment. 
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compromise between sensitivity and wire life. With higher operating 
temperatures more sensitivity is obtained, but there is a greater 
chance of contaminating or burning out the wire. So far, an overheat 
ratio of 1.4 has been used for all the experimental work. 

After adjusting the anemometer to the desired overheat ratio, 
data are taken for the pressure drop on the wind tunnel and the voltage 
output from the anemometer — in the range of velocities of interest. 
These data are then converted to velocity and voltage data, which in 
turn are used to determine the constants in the Col 1 is-Wi 1 1 i ams equation. 
This last step is done by a least-squares curve fit of the data to the 
equation, using a computer. Figure 18 gives the data, computer plot 
and resulting equation for the calibration of a wire. 

System Operation 

In order to take data, a through flow of air is first established 
in the test section, The velocity of this flow is adjusted to a level 
where the velocity is significant compared to the turbulent fluctuations, 
which is at a level of about 2 m/sec. In later work this velocity may be 
changed to establish its effect on the response. 

Next, the desired acoustic field is generated. This is done by 
installing the necessary acoustic drivers, tuning the General Radio 
1309-A function generator to a cosine wave at the resonant frequency and 
feeding this signal via a power amplifier to one driver, and, if two 
drivers are being used, adjusting the phase shifter to the angle required 
at the second driver. Then the power to the drivers is increased to a 
value in the vicinity of that necessary to generate acoustic pressures 
and velocities which are significant compared to the background noise 
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Figure 18. Hot Wi re Calibration Curve. 


and through flow velocity. This is at a sound pressure level of about 
160 dB. The phase shifter and the power to the drivers are then 
adjusted until the acoustic field has the proper characteristics. After 
that the desired distortion can be added by means of the Waveteck and 
Hewlett-Packard function generators. 

The proper characteristics of the acoustic field depend on the 
type field being studied. For example, to obtain a spinning transverse 
field with distortion, a spinning field without distortion is established 
first. The characteristics used to determine if the spinning transverse 
mode is present are as follows: 

•The amplitude of the sound is measured at numerous angular 
positions around the test section at a constant radial 
distance from the center. If the acoustic field is spinning, 
this amplitude should be constant. 

•The change in phase as the microphone probe is moved 
around the circle is also checked. If the acoustic 
field is spinning, the phase change should correspond 
to the amount of change in the angular position of 
the measurement. 

•The variation in the acoustic pressure along a diameter 

is measured. If it is spinning, (P ' /P ' ) should 

J max rms 

vary from one at the wall to near zero at the center 
with the decrease varying as a zero order Bessel 
function. Figure 19 compares such pressure data 
to the predicted Bessel function. 

• The change in phase along a diameter is measured. It 
should be constant from the wall to the vicinity of the 
center. In crossing the center, the phase should 
change by 180° and then be constant again to opposite. _ . 
wall. Figure 20 is a plot of the phase angle along a 
diameter for various radial positions. 




Figure 19. Sound Pressure Variation Along a Diameter. 
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Figure 20. Change in Phase Along a Diameter 



After the spinning transverse mode is established, the distortion 
is added by first locking the Wavetek and Hewlett-Packard function 
generators to twice the fundamental frequency, then gradually increasing 
their amplitudes and adjusting their phases relative to the fundamental 
signals until the desired pattern of distortion is achieved. This is 
checked by monitoring the anemometer and microphone signals on an 
oscilloscope. 

The theory correlates velocity and pressure at the same location, 
but experimentally both probes cannot occupy the same location. However, 
since the pressure varies in a predictable manner, as is illustrated 
in figure 19 , if the pressure at any location (except the center) is 
known, the pressure at any other point can be computed. Thus if the 
position of the hot wire is measured relative to the microphone and 
chamber, the microphone data can be modified to indicate the value of 
the pressure at the hot wire rather than the pressure at its actual 
location. 

The output of the anemometer goes to the Tektronix oscilloscope 
in order to remove E from the signal and to adjust the scale of the 
resulting signal so that it is within the 1-volt-rms limit on the input 
to the tape recorder. Also, the scale and range of the sound level 
meter are set so that the microphone signal is within the l-vo!t-rms 
limit of the recorder. When this has been done the two signals are 
recorded. A record is kept of the settings on all the instruments used, 
along with temperature and atmospheric pressure measurements, and 
calibration results. When this has been completed, the sound, microphone, 



and anemometer can be turned off, and the analysis can be done using the 
taped data. 

The analysis has been done in two ways, both of which have 
advantages; In the first method the time delay of the microphone signal 
is set manually at selected values and the corresponding correlation 
coefficients are all determined from the same section of the tape. This 
is accomplished by continually recycling the tape over the same section 
for each time delay. In this way a crosscorrelogram is obtained in 
which external effects are kept to a minimum. The room temperature, the 
barometric pressure, and the characteristics of the equipment are the 
same since the same data are analyzed over and over. However, this 
method is time consuming and tedious and results in data points rather 
than a continuous plot. 

The second method automates the procedure. A ramp function with 
a period of 2,000 seconds is used to drive the time delay unit. The 
amplitude of the ramp is adjusted so that it varies the time delay from 
zero to a value representing a complete cycle. The x-y plotter is then 
used to continuously plot a correlation coefficient versus the time 
delay. The disadvantage of this method is that the correlation coeffi- 
cient is being calculated at different sections on the tape and the data 
vary from one section to another. 

Results 

Data were recorded for several different types of sound field -- 
spinning with distortion, spinning without distortion, and standing 
transverse. Final results have not been obtained as yet, primarily due 
to delays in obtaining the proper amplifiers for the system. All data 
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taken with a spinning acoustic field contained distortions due to the 
irregularities in the sound field created by the temporary power 
ampl I f iers . 

In an attempt to account for the time delays introduced by the 
equipment, data were taken for the first standing tangential mode 
using a sine wave at the fundamental frequency. Pressure data were 
taken at the pressure ant inode and velocity data were taken at the 
center of the test section (the velocity antinode). In theory, these 
properties are related in a relatively simple manner. Figure 21 compares 
the theoretical and experimental results and indicates that there is 
either a 40° or a 220° shift in phase. The reason for the difference 
in amplitude of the two curves has not been accounted for. It may 
result from an error in the computations or may possibly be due to the 
method of dealing with E and E' experimentally. 

The correlation coefficient which is calculated experimentally is 


R = 


E' P' 

r p 
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The in-phase response factor can be obtained from R as follows: 
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Standing Tangential Mode . 
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where P is the barometric pressure and E is the time-mean of the voltage 
from the anemometer. 

Figures 22 and 23 show the experimental results obtained for the 
spinning transverse acoustic field with distortion as well as the 
analytically predicted results. As stated earlier, these results are 
not conclusive since they were obtained from sound fields distorted by 
the power amplifier. 
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Chapter 4 


CONCLUDING REMARKS 

The results obtained so far for both the analytical and 
experimental portions of the investigation Indicate that the idea of 
using a hot wire as an analog of vaporization limited combustion is 
feasible. The desired acoustic vibrations can be created and a response 
factor obtained. However, the instrumentation required In the exper- 
imental portion has been found to be much more complex than originally 
envisioned. The system as a whole is still in need of precise 
calibration. 

What has been done is only the p re 1 1 mi nary work needed to create 
a working system. To obtain information useful to design, a feedback 
system must be set up so the system responds naturally rather than 
having a predetermined pattern of vibrations forced upon it. To 
realize the full potential of the system, a method of time delaying a 
complex signal will be necessary. Also, the system needs to be enclosed. 
The sound levels created during operation preclude use during normal 
working hours, and those operating the system can only tolerate the 
noise for short periods, even when using proper ear protection devices. 
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